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Abstract: 

We discuss the high energy asymptotics in the next-to-leading (NLO) BFKL equa- 
tion. We find a general solution for Green functions and consider two properties 
of the NLO BFKL kernel: running QCD coupling and large NLO corrections 
to the conformal part of the kernel. Both of these effects lead to Regge-BFKL 
asymptotics only in the limited range of energy ( y = ln(s/ggo) < fas)" 1 ) and 
change the energy behaviour of the amplitude for higher values of energy. We 
confirm the oscillation in the total cross section found in Ref. in the NLO 
BFKL asymptotics, which shows that the NLO BFKL has a serious pathology. 
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1 Introduction 



The next-to-leading order corrections to the kernel of the BFKL equation [|TJ 
(NLO BFKL) have recently been presented in two papers @||, and have been 
formulated in a compact and elegant form by Fadin and Lipatov [[|. It turns 
out that the next-to-leading order corrections are large and they strongly modify 
the leading order result. The large value of the next-to- leading order corrections 
suggests that even high order corrections may be essential, leaving us without a 
reliable calculation of the parameters of the BFKL asymptotics except at very 
small values of the QCD coupling constant, which unfortunately are not presently 
attainable. 

Before arriving at a pessimistic conclusion, it is necessary to understand the 
qualitative alteration due to the next-to-leading order corrections that occurs in 
the BFKL asymptotics, i.e. we wish to know how the NLO BFKL Pomeron dif- 
fers from the LO BFKL one. Other questions which we would like to answer are: 
Does the NLO BFKL Pomeron still manifest itself as a Reggeon - like exchange 
at high energy as the LO BFKL Pomeron does? Can the NLO BFKL Pomeron 
be described as the diffusion cascading process in log kj_ where kj_ is the partonic 
transverse momentum? Can we calculate the NLO BFKL Pomeron in the frame- 
work of perturbative QCD or the nonperturbative correction would destroy our 
approach at so small value of energy that we cannot use the BFKKL asymptotic? 
In short, we would like to understand what the NLO BFKL Pomeron is ? 

This paper is an attempt to find a general solution to the NLO BFKL equa- 
tion, and to discuss the main properties of the high energy asymptotics that 
follows from this solution. The NLO BFKL kernel contains two parts: a con- 
formally invariant part and a running coupling part. In studying the running 
coupling effect we reproduce non-Regge type corrections to high energy asymp- 
totics suggested by Kovchegov and Mueller in Ref . [|J , using a method of solution 
proposed in Refs.|^|[^J. Ross found in Ref. 0] that the conformal part of the NLO 
BFKL kernel crucially changes the diffusion in log k± (kj_ is a parton transverse 
momentum ) which is a basic property of the LO BFKL Pomeron. We show how 
this affects the asymptotic behaviour of the scattering amplitude at high energy. 

The running QCD coupling in the BFKL equation has been studied for some- 
time, starting with the GLR paper ||, where the general solution to this problem 
was suggested. In Refs. §00 JUj gplpfg the different aspects of the 
problem were considered, but only now this can be done on solid basis of the 
NLO BFKL equation. The influence of the NLO on the BFKL diffusion was first 
pointed out and explored in Ref. f7|. 

Our paper deals with both problems on the same footing, and we give a simple 
and transparent derivation and discussion of the high energy asymptotics for the 
NLO BFKL equation. Our conclusion is rather pessimistic: we cannot avoid an 
oscillation in the total cross section at high energy in the NLO BFKL approach 
and, therefore, we do not think that the NLO BFKL equation can serve as a 
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basis for high energy phenomenology. However, we hope, that the experience 
with the NLO BFKL Pomeron will be useful for possible scenario of the high 
energy behaviour in QCD. 



2 A general approach to the NLO BFKL equa- 
tion 

The BFKL equation governs the high energy asymptotics of a single-scale hard 
process. As an example of such a process we can consider the scattering of two 
virtual photons ( with virtualities Q 2 and ) at high energy in the kinematic 
region where Q 2 — Qq ^> rn 2 ||14||, where m 2 a scale of the "soft" interaction. 
The cross section for this process can be written as 



a rr (s) = (2.1) 

where q and go are transverse momenta of gluons, &(Q 2 ,q 2 ) and ^(Q^q^) are 
impact functions that has been calculated in Ref.[T^]. It should be stressed that 
functions $ provide q 2 = Q 2 and q^ = Qq. G^q,^) obeys the BFKL equation 

uG u {?,$>) = 5®(q-q ) + J #(jK{U)G u tf,$ Q ) . (2.2) 

Kernel K(q,cf) can be written as a sum of LO and NLO kernels 

K{q^) = K LO (q, q) + K NLO (q,q>) . (2.3) 
This kernel has eigenfunction|| 

and corresponding eigenvalues 

d 2 q'K{q^)v f {q' 2 ) = u:{f) Vf {q 2 ) . (2.5) 



u)(f) has a form 

«(/) = (2.6) 
as(^)(x L °(f) + a s (fi 2 )x NL °(f)) ~ ^^- 1 bHq 2 / f , 2 ) X LO (f), 

7T 

where x LO (f) — 3-0(1) — 4>{\ + f) — "4>{\ — f), the explicit form of x NLO is 
written in Ref. 0, which has also symmetry under transform / — > — f as the LO 



part, and where b = jj^L_J^£ f or number of colours N c and number of flavours 
Nf. We use notation a s for as = 

One can see that the first term has a conformal symmetry while this symmetry 
is broken in the NLO due to the second term in Eq. (|2.6| ). fi 2 is the normalization 
point which is arbitrary in the NLO calculations. We take it to be equal to the 
value of the initial virtuality q 2 ( fi 2 = q$ ) without loosing any accuracy in the 
NLO approach. 

Following Refs. || we rewrite Eq. ( |2.6|) in the form: 



"(/) = 'fas{qi){x LU U) + as(q 2 )x NLU (f)) = j^Af) . (2.7) 

Here we define r = In — ~ \nas(q 2 ) and the running QCD coupling constant 
in leading log is equal to FI 

1 

or 

where A 2 is the position of the infrared Landau pole in running as- 

We want to stress that Eq. ( p.7|) coincides with Eq. ( |2.6|) in the NLO approx- 
imation. However, this form of the kernel is much more convenient in searching 
of a solution and we firmly believe it corresponds more to the general incorpo- 
ration of the effect of the running QCD coupling in the BFKL equation ( see a 
discussion of this point of view in Ref. || ). 

To find a solution to Eq. (|2.2|) we expand G^q 2 ) with respect to the complete 
set of eigenf unction of Eq. (|2.4|) , namely 

G^q 2 ) = (2.8) 

— g(w,f) (pf(q) = r == / — g{ujj)e\ 

ZITt Jq2 g2 Ja-ioc ZTTt 

where the contour of integration is situated to the right of all singularities of 
function g(u, f). 

Using Eq. (|2.5|) and Eq. (|2.7|) one obtains the following equation for function 

g(u,f) ■ 

- = r Q uj conf {f)g{ujJ) + roe-** . (2.9) 

The solution of homogeneous equation ( Eq. ( p.9|) without the last term ) can be 
easily found and it has the form ( see Refs. [Q @ for details ): 

g{u,f) = 9He"»^ m/(m ' • (2.10) 



1 We add factor ^\aas{q 2 ) for convenience but it does not affect the value of u){f) in the 
NLO approximation. 



4 



Function g(u) should be specified from initial or boundary conditions. The value 
of /o can be arbitrary since its redefinition is included in function g{uS). Unless 
it is specially stipulated fo = 0. 

We will show that the small values of / will be dominant in Eq. (|2.8|)in the 



wide range of large y = In t== . Therefore, we can expand ui C onf at small /, 
namely, 

UJconf = U L + Df + 0(f) , (2.11) 

where in the NLO ujl and D are equal to ( see Eq.(7) of Ref. @ ): 

u L = a s (q 2 ){ 2.772 - 18.3a 5 (g 2 ) } 5 (2-12) 
D = a s {ql) i 16.828 - 322a 5 (g 2 ) } • (2-13) 



One can see from Eq. (|2.12 ) how large and essential the NLO corrections 



are. They considerably diminish the value of lol which can be even negative for 
as(qo) > 0.152 and change the sign of D at as ~ 0.05. Note, that the positive 
value of D corresponds to diffusion in lng 2 . 

To solve this problem we need to find g{uj) in the general solution of Eq. (|2.10|) , 
which depends on the initial or boundary conditions. We find it very instructive 
to introduce two Green functions for the BFKL equation. 

1. 

The first one ( G r (y,r) ) satisfies the following boundary condition: 

G r (y,r): (2.14) 
G r (y,r = r ) = 5(y - y ) 

This Green function allows us to find us the solution of the BFKL equation for 
any boundary input distribution G in (y,q 2 = q 2 ) at q 2 = q 2 ( r = r ). Indeed, 
such a solution is equal to 

G(y,r) = [ dy G r (y,r)G in (y ,q 2 = q 2 ) . (2.15) 



Such a Green function is very useful for study of the boundary condition for 
the DGLAP evolution. Using G r (y,r) and Eq. (|2.15| ), we can investigate the 
^-dependence at q 2 = q 2 . We can distinguish two cases with different solutions: 

1. the integral over y depends mostly on properties of input function G in ; 

2. the integral over y is sensitive to the Green function. In this case we 
can claim that the energy behaviour of our boundary condition is defined by the 
BFKL dynamics. 

Therefore, this Green function ( G r (y,r) ) can provide us with an educated 
guess for the energy dependence of the boundary condition in the DGLAP evo- 
lution equations fL5|. 
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In the LO BFKL approach G r (y — yo, r, r ) is equal to[lj] 



G^°{y-y ,r,r ) 



Klr - r " ] . (2.16) 



We call this expression Regge-BFKL asymptotics since it has a power - like 
behaviour, similar to the exchange of the Reggeon. We would like to recall that 
the second factor in the exponent is different from a Reggeon contribution, and 
it has an origin in the diffusion in the log of transverse momentum of partons 
which is a typical feature of QCD. 

The question is how general is Eq. ( |2.16| ) and does it preserve the main char- 
acteristics like power behaviour and/or the diffusion in the log of transverse mo- 
mentum in the NLO approximation. 

2. 

The initial condition for the second Green function ( G y (y, r) ) can be written 
as follows: 

G y (y,r): (2.17) 
G y (y = y , r) = S(r - r ) 

It is obvious that one can find a general solution using this Green function 
if we have an input from experiment and/or nonperturbative QCD, namely, the 
dependence on transverse momentum at fixed value of y. Indeed, 

G(y,r) = J dr G y (y,r)G in (y = y ,r ) . (2.18) 

From Eq. ( [2.18| ) one can see that G y (y, r) defines the asymptotic of the scattering 
amplitude for a single-scale process. For example, this Green function gives the 
asymptotic behaviour of cx 7 * 7 * |L4 . 



The last remark in this general section: it is necessary to choose sufficiently 
large initial transverse momentum ( q% ) in order to safely apply the pQCD 
methods. It has been discussed in many papers H Hl6| |jl7H IfDfl |I^1 H that at high 



energies the BFKL diffusion in log of transverse momenta inevitabelly leads to 
the fact that small values of the transverse momenta become important (Bartel's 
cigar [[16]). It means that we cannot safely calculate the high energy asymptotics 
in the framework of pQCD. The criteria for being able to trust pQCD formulated 
in Ref.flTTfl (see also Ref. 0] for discussion in the case of running QCD coupling 
)is 

» £ uwh^ x sik (2 - 19) 

which suggests that we should take q$ as large as possible 



6 



In the LO BFKL approach this Green function has a similar form as G r , 
namely |3J 

G L y °(y - y , r, r Q ) = ^ £0 * _ ^ ^ } " . (2.20) 

One can see that both Green functions in the LO BFKL approach are similar 
and only differ in pre-exponential factors. 



3 G r (y,r) 

3.1 Solution 

This Green function has been calculated in Ref. but for completeness we will 
reproduce simple calculation to examine what happens to G r in the NLO. Sub- 
stituting Eq. (|2.11|) in Eq. fl2.10|) we find that the general solution of the homo- 
geneous equation is 



G{y,r) = — ^-g(u)^^)+tr (31) 

J a—ioo lit I Jfo-ioo LTX % 

The integration over / leads to Airy function [[J Ai f(^>)^ [r — ^7^0] V There- 
fore to satisfy the boundary condition of Eq. fl2.14|) we have to choose a function 
~g{u) = Ai' 1 ((^ [r - >„])• 

Finally 0], G r (y - y , r, r Q ) is equal to 



G r (y-yo,r,r ) = J- / — e»(v-v°> } r ° D \ ^4-. (3.2) 

3.2 Regge-BFKL asymptotics 

To recover the Regge-BFKL asymptotics we assume that the arguments of both 
Airy functions in Eq. ( |3.2| )are large: 

&Hr-^r ] = (^r [%- A] »i ; 

(*)Mro->o] = (*)ir [A] »1; (3.3) 



where 5r = r — r and A = uj — u^. 



, 3 

Using the asymptotics of Airy function |L9| v4i(z)| 2>0; | z |>i — > -^-e^^^ we 



obtain for positive D (D > 0): 



2z? 



G r (y-y ,r,r ) = / ^ e*^ , (3.4) 

J a—ioo Z7TI 
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where \I/ is equal to 



/ \ a / \ 2 fuJZ ( r Sr A 3 r A 3 \ 



We can calculate this integral using the saddle point method, in which the value 
of saddle point (A 5 ) is defined from the equation: 



dA 



|a=as 



. 



In our case this equation reads 

(y-vo) - 



1 ( ( Sr A* i r A* i 
Duj l \ r u L uj l 



We can easily find a solution to Eq. (|3.7|), assuming that 

Sr u L 



r A 



Eq. 



can be reduced to the form 



8r 



1 



1 e 

- K" 

4 



< 1 



1 

- ( K 



S \2 



2V-D (j/-y ) 

In the leading order with respect to k 5 , the solution to Eq 

Sr . c (Sr) 2 



2VD(y-y Q ) 



A s 



4D(y-y ) 2 



(3.6) 

(3.7) 

(3.8) 

(3.9) 
has the form: 
(3.10) 



However, we need to calculate the value of the saddle point with better accuracy 
to calculate the deviation from the Regge-BFKL behaviour. Considering Kq = 
4Dw L {y-y ) k e j n „ small, we find 

ro or ° ' 

Sr 



'A s 



2VD(y-y ) 



1 



(3.11) 



The integrand (see Eq. ( |3.5| ) ) for A = A s is equal to 



V(A s )=u; L (y-y )(l 



Sr 
2r~ n 



(Sr) 



+ 



1 Dul 



AD(y-y ) 12 r. 



(y-y ) 3 . (3.12) 



One can see that for positive D we have a good saddle point, and can evaluate 
the integral using the steepest descent method (SDM ). The result is 



1 



G r (y-yo,r) = — 



7r(Sr) 



qq \ 2D(y-y y 



(3.13) 
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The first three terms in Eq. ( |3.12| ) yield normal Regge-BFKL asymptotic form 
which is the same as in the LO BFKL approach. The difference is only in the 
second term which fixes the scale of the running QCD coupling in the expression 
for ujl (see Eq. (|2 . 1 2|) ), namely in the BFKL diffusion, the running QCD coupling 
enters at the scale q\ v = qqo). In other words we have to calculate u;£ un using 

= as(ql = qqo){2.772 - 18.3a s (g 2 )} (3.14) 



instead of Eq. fl2.12j ). This result was first obtained in Refs. HQ. We will show 
below that the same scale should be incorporated in the calculation of D. 

The last term in Eq. (|3.12| ) is the most interesting one, because it shows the 



violation of the Regge-BFKL asymptotics due to the running QCD coupling. 
It was suggested by Kovchegov and Mueller in Ref. || using quite a different 
method and, we hope, that our derivation presented here is more transparent. 
Eq. ( p.!2j ) shows explicitly that Regge-BFKL asymptotics is only valid in the 

limited range of y — y , namely, D( asuj L b) 2 (y — y ) 3 <C 1 or (y — y ) <C a s ' A 

EE- 

The question arises can we trust the term proportional to (y — y ) 3 - Indeed, 
we have made a lot of assumptions , in the derivation as well as using the steepest 
descent method (SDM). Therefore, we have to check whether all our assumptions 
are selfconsistent. They are: 



( W L 
\Dr t 



3 x A 1l° ^ ^ asymptotic of Airy function (3.15) 
k S — ^ 1 assumption used to obtain A s (3.16) 



JrSf Ia=as x [-^ < 1 self 'consistency of the SDM (3.17) 

r 2 

{y — Uo) < jn applicability of pQ CD (see Ref. [4]) (3.18) 



In addition, we have to check that ^||a=a s is positive, but it is easy to find out 
that it is the case for D > 0. However, it is an indication that D < it has to 
be considered separately. It turns out that the second equation of Eq. ( |3.15| ) is 
the most restrictive . Taking A s from Eq. ( |3.10| ) we see that it leads to 



» {^tbY x ^ (319) 




Substituting Eq. (|3.19|) in (y — y ) 3 -term in \1/ ( see Eq. (|3.12|) ) we obtain the 



estimate on maximum value of this term which we can guarantee in our approx- 
imation: 
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One can see that in wide region 5r > (as)~^, it is legitimate to keep this 
term and, therefore, we can conclude, that the running QCD coupling does not 
satisfy the Regge-type asymptotics. However, we have to understand this result 
better since, at first sight, zeros of the denominator in Eq. (|2.10| ) yield Regge 
asymptotics. 



3.3 High energy asymptotics ( a more detail analysis ) 

D > 

First, Airy functions have zeros only at the negative values of the argument, 



and their position can be found with good accuracy from the simple equation fl9 

, 3nn 37T , 2 
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(3.21) 



where z is the argument of the Airy function and n is arbitrary integral number 

(71 = 0,1,..."). 



Taking the argument of the Airy function in the denominator of Eq. (|3.2j) we 



obtain 0113 



/ D \ 3 „ 3™ 3tt , 2 . 

Here,A n = uj n — ui and one can see that all poles in u> are located to the 

2 

left from ujl ( u n < ujl and A n oc agOJL )• Therefore, we can legitimately 
calculate A n in framework of our approach. The whole structure of singularities 
in the ui - plane is as follows || |13| : 

1. the rightmost pole ojq is located to the left of ujl but, theoretically, very 
close to it. Namely, its position is 

( D \ 5 . 3tt 2 
uj = u L - u L [ 2 x ( — )3 . (3.23) 



2. for large n —>■ 00 



D \ 3 



Wnk»i -> ^1 x ^— ; (3.24) 

3. there fore, the solution has the infinite number of poles in tu-plane to the 
left of ujq , which accumulate at uj = at n — > 00. 

This picture of singularities, justifies our saddle point calculation, since the 
position of the saddle point turns out to be shifted to the right of u>l- Note, that 
the contour of (^-integration is chosen to be located to the right of all singularities 
of the integrand. 
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We can evaluate the integral in a different way, namely, calculating each pole 
separately. In this case we have the asymptotic form 



G r (y-y ,r,r ) = £ e^^~^V(u; n ,r) , (3.25) 

n=0 



where V(u n ,r) is the u n - pole. This series was investigated in Ref.[]T3|] where 
it was shown that the saddle point approximation effectively describes the sum 
over n at sufficiently large n. The right most singularity ojq gives a suppressed 
contribution because of smallness of its residue, while uj n — > at large n, but 
the residue is rather big. 
D < 

In this case the general solution has the same form of Eq. ( |3.2|) but the struc- 
ture of singularities changes crucially: 

1. all poles u n are located to the right of ul, 

2. the position of the leftmost pole is 

( D \* ,371-, 2 . 

u = uj l + u L \ — 2 x (^) 3 ; ( 3 - 26 ) 

3. at large n u n — > oo accordingly the following expression: 

\D\ /37rn\ 2 

w„| n >i = V x hr • ^ 3 - 27 ) 

r \ z / 

In such a situation we cannot use the saddle point method and should rather 
analyze Eq. ( |3.25| ) with u n given by Eq. ( |3.27| ). It is easy to see that Eq. ( |3.25| ) 
can be reduced to the form: 

G r (y-y ,r,r ) = ^ eT 2 r ° . (3.28) 

n>l 

To evaluate this sum looks hopeless, at least, we do not see how it is possible 
to do . However, we have to go back to our derivation of solution (see Eq. ( |3.1| )), 
namely, to the integration over / which led to the Airy function. The value of the 
typical / ( f s ) in this integral can be evaluated using the saddle point approach 
and it is equal to 

f s - ^W\ {r - r ^ ] ' (3 ' 29) 

which is of the order of n for u n . Therefore, for large n we cannot trust the 
solution of Eq. ( |3.2|) and have to generalize the solution including the next term 
in expansion of u ccm f(f) in Eq. (|2.11|) , namely, 

Uconf = ujl + Df 2 - Bf\ (3.30) 
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where B was calculated in Ref.]7| 

B = -a s (q 2 o) {64.294 - a 5 (g 2 ) 2756 } . 



(3.31) 



One can see that B > everywhere, except of the region of very small as- 
One can see that Eq. ( |3.30| ) can be written as 



OJ. 



conf 



*L ~ fi(/ 2 "/o 2 ) 2 



where 
and 



wl = ujl + B f 



(3.32) 



f 2 _ _ \D\ . 
~~ 2B ' 
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Introducing a new variable f = fo + v and integrating in Eq. fl2.10|) from 
fo we obtain 

G(y,r) = (3.33) 



a -ioo 2ixi Jf-ioo 2ixi y ' [ V 2B 1 



Integration over v in Eq. ( |3.33|) gives the Airy function Ai ( (^757) 3 [r — ^7 ?"o] 



with 

£>' = 2\D\ . (3.34) 
Finally, the solution is very similar to Eq. ( |3.2| ) with new D' of Eq. ( |3.34| ), namely 

Gr(y - Vo,r,r ) = (3.35) 



, r^du \D\ } 
r Ja-ioc 2m \ V 2B 



Therefore, Eq. (|3.12|) gives the answer in this case, for the saddle point value of 
\]/ with D = D' . The Green function G r (y — yo, r, r ) has the form: 



G r {y - yo,r,r ) = * 



2D'{y-y Q f [ V 2B 1 



The unpleasant fact which is a direct consequence of the NLO BFKL is the 
oscillation, which has first been found in Ref. [0. This is a serious defect of the 
NLO BFKL, since G r is proportional to the total cross section. Certainly, we 
cannot have faith in such an approach which leads to negative total cross section. 
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3.4 Asymptotics at ultra high energies 

We can calculate the asymptotics of this Green function ( G r ) even at very 



high values of energy Indeed, Eq. (|3.18| ) which restricts the value of y — yo does 



not affect G r since the nonperturbative corrections have been absorbed in Gi n in 
Eq. (|2.15|) . Eq. C P .15 ) gives a nice example how we can separate the unknown 



nonperturbative input distribution ( G in ) and perturbative BFKL Green func- 
tion. Unfortunately, we are not able to do so in the case of the second Green 
function ( G y ) for which Eq. ( |3.18j ) is a restriction ( see Ref. M for details ). 

_ 5 

It is easy to see that at very large value of y — yo ^> a s 3 the inequalities 
of Eq. ( |3.15|) - Eq. ( |3.18|) are violated and actually we cannot use the asymp- 



totic expression for the Airy functions both in numerator and in denominator. 
Therefore, the only source of the asymptotics is the zero of the denominator at 
uj = u Q , where (see Eq. (|3.23|) ) 



= ^ {tY • (3 - 37) 

Closing contour in u - integration on this pole we obtain 

G r (y - Vo,r,r ) = (3.38) 

" TAi({^-)Ur - A-A x joiv-vo). 

s V t o L) 8 / 



8D 2 



Therefore, we obtain the typical Regge asymptotics without any diffusion in the 
log of transverse momentum. The first factor is the Reggeon residue while the 
second one is the Reggeon propagator. 

For the case D < one can see that the asymptotics will be given by 
Eq. ( |3.37| ) with obvious substitute u L — > uj (see Eq. ( |3.32| ) ) and D — > D' = 
2\D\ ( see Eq. (ggg ) ). 

One can see, that Regge behaviour depends only on initial virtualities Qo{r ) 
and corresponds to the pole in the angular momentum with the position which 
is independent from any characteristics of the "hard" processes. We hope, that 
this example can be instructive for high energy phenomenology. At least, it gives 
an answer to the question: why and how the Regge-like behaviour could appear 
for the initial condition for the DGLAP evolution equations. Recall, that such a 
behaviour is heavily used in the solution of the DGLAP evolution equations. 

4 G y (y,r 
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4.1 A general solution 

To find G y (y, r) we return to a general equation of Eq. (|2.9| ). The inhomogeneous 
term in it corresponds to Eq. (|2.2|) or, in other words, the inhomogeneous Eq. (|2.9| ) 
is written for G y (y,r). One can easily find the solution to this equation, taking 
the solution of the homogeneous equation ( see Eq. (|2.10| ) ) but considering g as 



a function of to and /. Substituting it back to Eq. ( |2.9| ) we obtain the following 
equation for g(u, f): 

dg{u, f) _ /ro a f!)u> conf (f')df' v 

— to — — — = r e 3 °e" Jf o . (4.1) 

°y 

Finally, a general solution can be reduced to the form: 

d w rf+ioc d f rQ roc d f> 

Gy{y -y ,r,r ) = — / _ — / — (4.2) 



ioo 2 7r z Jf-ioo 2ni u Jf 2iri 

e u(v-Vo) + fr - f'ro-% ff Uconf(f")df" 

The integral over u can be evaluated and we obtain the general solution in the 
form that we will use in this paper; 

G y (y-y ,r,r ) = r [ [ %^6{f - f)^~^ (4.3) 

J J 2m 



I J 2 J(y-y )r f f U conf (f")df" 



f 

It is easy to check that G y (y — yo, r, r ), given by Eq. ( |4.3|) , approaches 5(r — r Q ) 
at y -> y . 

Since at y — yo ^> 1 the argument of Iq is large, we can use the asymptotics 
of the modified Bessel function, namely, Iq(z) | 2 >i — > ^T^e 2 - To get a more 
compact answer we introduce: (i) new variables /~ = /' — / and f + = f + f] 
and (ii) the integral representation for 9 - function 

O(f-) = / ^ e . (4.4) 



The final result is 



G,(y-y ,r,r ) =r / / / |^ (4.5) 



/+ + /- 



e 
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4.2 Regge-BFKL asymptotics 



We assume that both / and / + are small enough to use Eq. (|2.11|) for uj con f(f). 
We can see that the answer has a form: 



G y (y - Vo,r, r ) 



df + df d\i y 



(y,r,f+,f-) 



(2ttz) 2 /i 

and taking explicit integration over /", * can be written as 

2^ L r (y-y ) y/J= ■ {1 + £- L [\U + ? + lit?]}- 



(4.6) 



(4.7) 



We evaluate the integrals over / and / + , using the steepest descent method in 
which we have the following equations for the saddle points : 



\f-=f- 







<9* 

df+ 



1/+=/+" 



. 



(4i 



Eq. (|J) gives 
[i + 



r + rp 
2 



^to-rt x {1 + £-[(f + ' S ) 2 + Ur S ) 2 }} = 0(4-9) 



VP 



r - r 
2 



8oj l L \J I 1 3 

Ju L r (y - y )VT X iflT f + ' S = 



The leading order solutions of the above equations are 

2 



/o 

Jo 



.5 



2^; \/^ r ° (f-w) Vfo 

Keeping the next order term we have 



(4.10) 

(4.11) 
(4.12) 



f+ ,s = 



y/u L ro(y-yo) \ ^ f1 _D_\(f+,S\2 



r + rn 
2 



x {i + ife[(/o + ' a ) 2 + iUT 5 ) 2 ]}; (4-13) 

— {1 - 



y/u L r {y-yo) Vfo' 



£-J(fo' S ) 2 + !(/o^) 2 ]}-(4.14) 



Substituting Eq. (|4.13|) and Eq. Q4.14|) in Eq. (|4.7| ), we obtain for the saddle value 
of* 



^ (v-vo. 



{Srf 



J_ n r+rp 
12 2r 



(^^) 2 (y-Z/o) 3 - 



+ 



(4.15) 
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Here, we have taken the integration over //. Indeed, it is easy o see that the 
saddle point value of (^y,r,f +,s ,f~' s ^j — > — oo when // — > — oo. Therefore, 
we can close contour of integration over and take a residue at /i = 0. 

In Eq. ( |4.15|) one can see the power increase (Reggeon - like) with energy ( 
the first term ) and BFKL diffusion in log of transverse momenta ( the second 
term). It is interesting to note that in both terms as well as in the third one, the 
running QCD coupling constant depends on the scale q 2 ev = qqo(r ev = rs 7± r ). 
Namely, coupling constant at this scale enters in calculation of both ujl and D. 

For both Green functions we obtain a term proportional to (y — y^) 3 which 
was first found in Ref. [g] Q. 



4.3 Regge-BFKL asymptotics for D < 

One can see from Eq. ( [4. 1 1[ ) that the saddle point value of /" does not depend 
on the value of D. It suggests the following logic of approach: we expand \I/ in 
a general solution of Eq. ( |4.5| ) in the region of small /~ while keeping / + to be 
arbitrary large. In such an approximation we have 



-f- ^conf (f)df' 1 



f + +r\ (f + -r 



df- 2 



- { ^conf ~ + U conf } 



"conf (£) + I J' con S {^\ (D' , (4-16) 



wVlOT-o I }" — rf2a; eonj(/) 

wnere to con f — — -^j — • 



The value of / in the saddle point is equal to 



f-rS = (4.17) 



, f + \ , , I. 5 u/'cnf (g) 

r uJconf — (y - y ) 1 + ^ 77+v ( fo 



,S\2 



where 



fo' S 



r + r\\ roUJconf ^ ~ ^ ' ^ 4 ' 18 ^ 

\I> for this saddle point is equal to 

y(y-y ,r,f+,r> s ) = (4.19) 



2 The numerical coefficient in front of this term (^) is different from the coefficient calculated 
in Ref. (±) due to different definition of D ( D(our) = 4 D(RefM) ). 
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Sr 2r (f+ \ 1 t/W (4) , ,- 



,S\2 



The position of the saddle point for / + integration is located near to the minimum 
of uJcmfi^)- If we denote the position of minimum / = Ref + i/m/o we can 
obtain the result after taking the integral over / + using the saddle point method. 
In practice as we have mentioned in the NLO BFKL the minimum occurs at 



fo = ±i V§ ( see Eq. flOg) ) for D < and f = for positive D. 



2B 

The result is 



UJ, 



conf 



(*) 



«*^" cm/ (f) (^) 

where 



G y (y-y ,r,r ) = (4.20) 
- x cos{2Imf 5r) x e *(f-fo,r-,r ,/+^/"- s ) ^ 



/ i c o\ ^conf (¥) 2r 
* (y - 2/o, r, r , /+' s , r < 5 ) = ^ 0/ - l/o 4.21 

(5rf 1 „ ^ ^on/(f) , 
+ ^ con/ I y ) , r „ +M2 (y - Jfo) • 



' r +T" \2 



2u/' C(W /(f) ( y - y 

5 Summary 

In this paper we analyze the prediction for high energy asymptotics that emerges 
from the NLO BFKL equation. Below we summarize the results of our analysis. 

1. We found two Green functions ( see Eq. ( |3.2p and Eq. ( |4.5| ) )that govern 
the asymptotic behaviour of the scattering amplitudes at high energy with two 
different initial conditions: at fixed virtuality Q 2 = Ql and at fixed energy (x). 

2. For sufficiently small values of energy ||[|J 



y = ln(s/qq ) < a ' 



5 

S 



in both Green function we found the Regge - BFKL asymptotics 

G oc e "^ y - yo) - ^v- 2 vo) (5.1) 

However, parameters u a and D a in Eq. fl5.1|) turns out to be different for dif- 
ferent sign of D in the NLO expression for the BFKL kernel (see Eq. ( |2.11| )) 

D > D < 

Lu a = ^- o u L (seeEq. u a = ^ [u L + g] (see Eq. Q£3D) 

D a = -^n- D D a = 2-2*>-\D\ 

r+ro r+ro i i 
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3. We confirm the result of Ref. |3J that the first corrections to the Regge- 
BFKL asymptotics lead to an additional term in the exponent of Eq. (|5.1|) , which 
has a form: 

(5 - 2) 

We showed that it is legitimate to consider such a term in the framework of our 
approach. 

4. Unfortunately, we found the oscillation in r — r for the asymptotics in 
the NLO BFKL in the huge region of initial virtualities 5o( r o) where the NLO 
generates D < ( see Eq. flOD ). This oscillat ion was first suggested in Ref. 



|. Indeed, the asymptotics has a factor ( see Eq. (|4.20| ) 




G y oc cos \ 2(r-r )\^ ■ (5.3) 



We consider this result as a challenge for the experts in QCD since it is a strong 
indication that the NLO BFKL approach has a pathology. We do not think that 
an additional integration in Eq.(2.1) will rule out this oscillations.lt should be 
stressed that the exchange of the BFKL Pomeron has a clear physical meaning 
as the asymptotics of the colour dipole - dipole scattering ( see Ref.|20| ). It is 



difficult to believe that anyone would cope with a negative total cross section for 
such a process []. 

5. Our solution of Eq. ( ^2|) suggests the asymptotic behaviour in the region 
of large values of Q 2 . Indeed, the Regge-BFKL asymptotics as well as the term 
which violates it come from the saddle point approximation in the region of 
uj « u a , but we find that uo a falls at r ^> 1, while the contribution from 
the zeros of the denominator in Eq. (|3.2f ) do not depend on virtuality Q 2 {r). 
Therefore, at large Q 2 the rightmost singularity in the u integration will be the 
position of the zero 

*o = o, a - (^) (^) 3 • (5.4) 

This singularity leads to Regge - behaviour 

G r oc e Wo(s/ - w,) 

in spite of the fact that in the huge region of intermediate values of Q 2 we have 
a much more complicated behaviour. 



3 It should be stressed that it was first demonstrated in Ref.[£lJ that the NLO BFKL 
corrections to the anomalous dimension lead to the negative probability at exceedingly small x. 
Perhaps, this result is correlated with ours, but the difference is that we observe our oscillations 
in the region of applicability of the BFKL approach, and, therefore, we cannot follow the "wise" 
advise of Dokshitscr |22]:"Let DIS structure functions in peace!" 
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This observation can be used as a justification of the Reggeon-like behaviour 
of the structure functions in the region of small x, used for the initial condition 
for the evolution equations. 

6. We cannot trust the perturbative QCD calculation for the Green function 

2 

G y since at y — y < tjj ( see Ref. [|J ) the contamination from the nonpertur- 
bative QCD region is rather big. We think, this is a principal difference between 
G r and G y , that in G r all nonperturbative corrections can be included in the 
input distribution Gi n in Eq. ( |2.15| ) while G r can be calculated in the framework 
of perturbative QCD. 
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